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Abstract. Spectral analyses of extreme horizontal branch (EHB) stars indicate a 
trend to supersolar heavy-metal abundances that is increasing with effective tempera- 
ture. The sdOB EC 1 1481-2303 is such a star, located at the very hot end of the EHB. 
In contrast to other sdOB stars, it is exhibiting an extraordinary flat UV continuum. 
This may be the result of extreme metal-line blanketing. We present preliminary results 
of a spectral analysis by means of static as well as stratified NLTE model atmospheres. 



1. Introduction 

Subdwarf B (sdB) stars form a homogeneous class of extreme horizontal-branch (EHB) 
stars (20 000 K 5 Teff ? 40 000 K, 5 < log^ < 6) that have evolved directly from the red 
giant branch. They are core helium-burning stars with hydrogen envelopes thin enough 
to prevent hydrogen burning. Their spectra show broad B aimer lines up to a principle 
quantum number of « « 12 and mostlj no helium lines. sdOB stars show similar spectra 
but additionally a weak He u A 4686 A line. They are found at higher Tes (> 40 000 K). 

In the last few years, many abundance analyses of subdwarf stars were performed. 
E.g. Edelmann (2003) analyzed a sample of 146 sdB and sdOB stars and found con- 
stant abundance values for similar T^s and log g but a general trend of an increasing 
helium abundance with rising T^s. Three stars showed extremely overabundant iron- 
group elements (up to 32 000 x solar). Geier et al. (2010) performed an abundance 
analysis of 139 sdB stars and detected a trend of increasing heavy-element abundances 
(generally supersolar) with Jeff whereas iron itself is not affected. This analysis showed 
that the super metal-rich stars of Edelmann (2003) are not that peculiar considering the 
supersolarity of heavy-metal abundances. In fact such values can be explained with 
atomic diffusion (equilibrium between gravitational and radiative forces), that becomes 
more effective at higher and confirms the trend found by Geier et al. (2010). 

There are still some features that cannot be explained. If pulsations are consid- 
ered, iron is believed to play a role in the driving mechanism. Therefore precise deter- 
mination of iron-group abundances are important to understand the difference between 
pulsating and non-pulsating stars and to explain the pulsation processes. 

One very peculiar star is EC 11481-2303, a relatively hot sdOB star {Tf.f\- = 
55 000 K) that shows extreme overabundances of iron and nickel and exhibits a pe- 
culiar, flat UV continuum that is not explained. Its atmosphere seems to be dominated 
by iron-group elements and therefore, it is an ideal candidate to test static and stratified 
models as well as stellar evolutionary theory. Spectral analysis of such a metal-rich star 
is still very difficult, mainly because of the lack of reliable atomic data. We summarize 
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previous analyses in Sect. 2. Our own analysis procedure is described in Sect. 3 consid- 
ering static as well as stratified model atmospheres, and the interstellar Une-absorption 
spectrum. Ways to access the TMAP model atmosphere code and its atomic data are 
shown in Sect. 4. We conclude with our preUminary results in Sect. 5. 

2. Discovery and Previous Analyses 

Discovered in the Edinburgh-Cape Blue Objects survey, EC 11481-2303 is a relatively 
bright source with my = 1 1.76 (Kilkenny et al. 1997) and was classified as a DA white 
dwarf with a companion 6'.'6 away. 

Stys et al. (2000) performed a first spectral analysis with LTE models considering 
hydrogen and helium only. They used optical SAAO (Aug 15, 1995, 1.9m telescope, 
3350-5450 A, resolution 3.5 A) and lUE spectra (Jul 14, 1993, 1150-1978 A, reso- 
lution 7 A for SWP48111 and 0.1 A for SWP48112). From the optical spectra they 
determined Teff = 41 790 K, log g = 5.84 and He/H=0.014 (number fraction), but were 
not able to fit the UV continuum with these models. From the UV spectra, they deter- 
mined lognm = 20.5 and concluded that EC 11481-2303 must have the same distance 
(500 pc) and reddening (0.06) as the nearby (distance 4?4) star HD 104377 for which 
the same log rim and similar line-of-sight elements and their velocities were determined 
by van Steenberg & ShuU (1988). Considering this distance, Stys et al. (2000) deter- 
mined an absolute visual magnitude of My = 3.2 for EC 1 1481-2303 which is typical 
for a subdwarf star. Therefore EC 11481-2303 was re-classified as subdwarf star. 

In order to explain the peculiar flat UV continuum, Stys et al. (2000) consid- 
ered the analysis of the DAB white dwarf GD 323 performed by Koester et al. (1994). 
GD 323 shows a similar flat continuum-flux shape. They considered models with ho- 
mogeneously mixed composition, chemically stratified constitution, a He spot at the 
surface, and a composite of a DA and a DB white dwarf. Only stratified models were 
able to fit the flux shape, but those models did not match the B aimer lines. Therefore the 
flux shape was still unexplained and Stys et al. (2000) concluded that EC 11481-2303 
is worth a follow-up study. 

Within the framework of the ESO Supernovae Type la Progenitor Survey (SPY) 
project (Napiwotzki et al. 2003), a sdO sample was observed with the UVES spectro- 
graph at the VLT and analyzed by Stroeer et al. (2007). EC 11481-2303 was excluded 
from this analysis among other things because of its UV flux shape. However they 
suggested that EC 1 1481-2303 has a higher Teff and lower helium abundance than de- 
termined by Stys et al. (2000). This SPY spectrum was used also by Jordan et al. (2007) 
to investigate magnetic fields in DA white dwarfs, but no hints of a magnetic field were 
found in the observations of EC 11481-2303. 

Rauch et al. (2010) used this spectrum together with the previously mentioned lUE 
and additional FUSE observations. With models calculated using the Tubingen NLTE 
Model-Atmosphere Package (TMAP, Werner et al. 2003, Rauch & Deetjen 2003) con- 
taining the elements H, He, C, N, O, and (Ca-Ni), they determined T^g = 55 000 K, 
\ogg = 5.8, He/H=0.0025 by number, and upper limits [C] = -0.875, [N] = -0.875, 
and [O] = -0.875 from the optical spectrum ([X] denotes log [abundance / solar abun- 
dance] of element X). To match the UV flux shape, they concluded that the abundances 
for the iron-group elements must be at least ten times solar and that nickel seems to 
be even more overabundant. They suggested that diffusion calculations should be per- 
formed in a next step and the iron-group elements should be considered individually 
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(they used a generic model atom). The reddening should be redetermined, and inter- 
stellar lines should be modeled. All this was picked up for this analysis. The exact 
procedure is described in the next section. 

3. Analysis procedure 

In the following, we introduce the model-atmosphere code for static model atmo- 
spheres, the statistical treatment of the iron-group elements, and the program for strat- 
ified atmospheres. Then we describe the Une-fitting procedure and the program for the 
interstellar lines. 

3.1. TMAP - The Model-Atmosphere Code 

TMAP^ has been developed since the 1980s and is well tested in many spectral analyses 
(e.g. Wassermann et al. 2010). It assumes hydrostatic and radiative equiUbrium and 
can consider plane-parallel or spherical geometry. TMAP provides reliable results for 
temperatures between 20 000 K < T^q < 200 000 K and 4 5 logg 5 9 and up to 1500 
NLTE levels can be included. An input for this code is model atoms from the Tubingen 
Model Atom Database (TMAD^). This database was built up together with the model 
atmosphere code, therefore its format is TMAP compliant. It is continuously updated. 
Presently it contains data for the elements H, He, C, N, O, F, Ne, Na, Mg, Si, S, Ar, 
and Ca. Data for other elements are in progress. Model atoms for model-atmosphere 
calculations as well as for SED calculations (that include fine-structure spUtting) can 
be downloaded. 

3.2. IrOnIc - Statistical Treatment For Iron-Group Elements 

Iron-group elements (Ca - Ni) consist of hundreds of milUons of lines due to their 
electron configuration (partly filled 3d and 4s shells). Therefore a statistical treatment 
is necessary to consider all lines but to get a reasonable number of NLTE levels. This 
treatment is implemented in the Iron Opacity and Interface code (IrOnIc, Rauch & 
Deetjen 2003) that divides the excitation-energy range of an ion into several bands and 
combines all atomic levels within one band to superlevels. The line-absorption cross- 
sections are sampled to superlines, ranging from one band to another or within one 
band. With this treatment all levels are considered, but the number of NLTE levels is 
strongly reduced. 

The program uses POS and LIN lines from Kurucz' line lists^ as input. POS lines 
contain only observed (good) wavelengths, LIN lines include theoretical as well as 
observed wavelengths and there are therefore about 300 times more of them. 

At first, we calculated a HHeCNO model atmosphere with the values determined 
by Rauch et al. (2010). To re-determine the reddening, we used the FUSE and lUE spec- 
tra as well as optical and infrared magnitudes from SIMBAD"^ (Table 1). We normalized 



'http;//astro.uni-tuebingen.de/~TMAP/TMAP.html 
^http://astro.uni-tuebingen.de/~TMAD/TMAD.html 
^ http ://kurucz. harvard, edu/ atoms.html 
'*http://simbad.u-strasbg.fr 
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Table 1. App arent magnitudes of EC 1 148 1 -2303 from SIMBA D. 

B V J H K 

11.49 11.76 12.512 12.669 12.765 
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Figure 1. The UV spectrum of EC 11481-2303 compared to a model contain- 
ing HHeCNO (blue, thin line) and hydrogen and iron-group elements with 10 x so- 
lar iron and 1000 x solar nickel value (red, thick line). Gaussians of lA and 3 A 
(FWHM) are applied to the FUSE observation and synthetic spectrum for clarity. 
The synthetic spectra are normalized to the reddest magnitude (Table 1). A redden- 
ing of Eb-v = 0.03 (0.01 for the HHeCNO model) is applied. 

the synthetic flux to the magnitude and applied the law of Fitzpatrick (1999). A 
value of Eb-y = 0.06 as determined by Stys et al. (2000) is too high. With our models 
Eb-w ^ 0.02 ± 0.01 fits best. 

For this analysis, we started with a model atom including hydrogen and all iron- 
group elements individually. We adopted T^s and logg from Rauch et al. (2010) and 
increased the abundances of single iron-group elements up to 1000 x solar. With these 
models, we compared our synthetic spectra with the observations (FUSE: B0540901 
from May 2001 and the previously mentioned lUE and UVES observations) and found 
that only the elements Ni and Fe are able to decrease the flux shape notable in the 
relevant regions (Fig. 1). We tested combinations of these elements with different su- 
persolar abundances. The flux is reproduced best by a model with 1000 x solar nickel 
abundance and 10-100 x solar iron abundance. The increase of the other iron-group 
elements to supersolar values to this model did not reduce the flux significantly. 

3.3. NGRT - The Interplay of Radiative Levitation and Gravitational Settling 

The atmospheres of sdOB stars are dominated by diffusion. We used the Next Gen- 
eration Radiative Transport (NGRT, Dreizler 1999) code that considers this effect. 
The theory was originally developed by Chayer et al. (1995) and modified by Drei- 
zler (1999). NGRT uses the TMAP model as input. In an iteration cycle, the radiative 
acceleration is calculated in a fixed atmosphere structure, then the radiation field and 
the atmospheric structure are computed again, etc. 
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Figure 2. Same as Fig. 1. The model considers diffusion and E^-v - 0.025. 

With the model containing hydrogen and the iron-group elements with solar abun- 
dances as start model, the NGRT calculation provided the abundance distribution shown 
in Fig. 2. This model reproduces the observation in many parts of the wavelength range. 
It also confirms the results determined with the TMAP model, both yield (extreme) su- 
persolar abundances for iron and nickel (Fig. 3) and show a flat UV continuum. 

3.4. Detailed Modeling of the Line Profiles 

The FUSE and lUE observations show many deep and broad absorption features that 
seem to stem from blended finer lines. With a pure HHeCNO model, only a few lines 
can be identified. Most lines seem to originate from iron-group elements. 

In order to model the line profiles, we compared the best-fit static and stratified 
models with the FUSE and lUE observations. Surprisingly, models and observations 
seem to have nothing in common although some individual features match the position 
of the observed lines considering the TMAP model (Fig. 4). The same is valid for 
individual regions of the spectrum. E.g. from 996 - 1004 A the model can reproduce 
the observation while nearby regions of the spectrum do not match at all (Fig. 4). 

The theoretical lines are much narrower than observed. This could be a hint for 
stellar rotation. The rotational broadening smears out narrow lines and they appear as 
one broad line. In the analysis of Rauch et al. (2010), the model shows a slightly higher 
emission reversal of the Hor line than the optical observation. With a rotational veloc- 
ity of Vrot = 30 km/ sec applied to their model (Fig. 5), this emission feature matches 
the size of the observed line. With this Vrot, more features in the FUSE observation 
can be reproduced. Rauch et al. (2010) included iron-group elements with a generic 
model atom into their calculations. In the fit to the Balmer lines, their model shows the 
so-called Balmer-line problem Napiwotzki & Rauch (1994). Thus, the heavy-metal 
abundances are still too low or not enough metal opacity is considered. 

We calculated models with individually increased abundances of the iron-group 
elements. The comparison with the observation shows that no line can be unambigu- 
ously identified. The lines mainly stem from Fe and Ni, but cannot be reproduced 
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Figure 3. Depth-dependent abundances of the single iron-group elements in our 
NGRT diffusion model. [X] denotes log [abundance / solar abundance] of element 
X. The solar abundance is displayed by the dashed line. 



perfectly. The reason is most likely that the Kurucz' LIN lines bear some uncertain- 
ties. We use Kurucz' line lists as input and only a few wavelengths are observationally 
confirmed (POS) for the high ionization stages that we have to consider in the case of 
EC 11481-2303. About 300 times more lines are only theoretically calculated (LIN) 
with uncertainties of wlO% in wavelength and wl5% in f-values. All LIN lines have to 
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Figure 4. Section of the FUSE observation compared with the best-fit TMAP 
model. A Gaussian of 0.06 A (FWHM), a radial velocity of -17.5 km/ sec, and a 
rotational velocity of 30 km/ sec are applied to the synthetic spectrum. A neutral 
hydrogen column density of 3 xlO^"/ cm^ and a reddening of Eb-v = 0.02 are con- 
sidered. 
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Figure 5. The theoretical Balmer-line profiles with Vrot = and 30km/ sec (green, 
thin and red, thick, respectively) compared with the observation. Additionally our 
best-fit model with 10 x solar Fe and 1000 x solar Ni abundance is displayed with 
Vrot - 30 km/ sec (blue). 



be included in the calculations because otherwise a huge amount of opacity is missing. 
But this could be the reason that the lines do not fit perfectly or seem to be shifted a bit. 



3.5. OWENS - Modeling Interstellar Absorption Lines 

To model the interstellar line-absorption spectrum, we use the program OWENS. It 
can model different interstellar clouds with different temperatures, radial and turbulent 
velocities, and column densities. The resulting ISM spectrum is combined with our 
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TMAP spectrum and some of the strong lines can be identified as interstellar (Fig. 6). 
But still the majority of the lines is of photospheric origin. 
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Figure 6. Section of the FUSE observation compared with the best-fit TMAP 
model (red) and the normaUzed OWENS ISM model (blue). The TMAP model has 
the same applications as in Fig. 4. 



4. TheoSSA - Virtual Observatory 

The German Astrophysical Virtual Observatory (GAVO^) aims to make astronomical 
data of all kinds accessible worldwide. It is part of the International Virtual Obser- 
vatory Alliance (IVOA^) and is a contact point for everybody needing assistance with 
publishing his or her own data or using tools or services. 

Within the framework of a GAVO project, TMAP was made accessible: 

1) Already calculated spectral energy distributions (SEDs) can be downloaded via the 
service TheoSSA (http://dc.g-vo.org/theossa). 

2) Individual models can be calculated with the web-interface of TMAP with the service 
TMAW (http://astro.uni-tuebingen.de/~TMAW/TMAW.shtml). 

3) Ready-to-use model atoms can be downloaded and used for every code or tailored 
individually from the Tubingen Model Atom Database 

{TMAD, http://astro.uni-tuebingen.de/~TMAD/TMAD.html). 

We will include the final SED calculated for this analysis into TheoSSA. 



^http://www.g-vo.org/ 
'http://www.ivoa.net/ 
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5. Preliminary Results and Conclusion 

We calculated static and stratified model atmospheres to reproduce the peculiar, flat 
UV continuum flux shape of EC 11481-2303. The elements H, Ca, Sc, V, Cr, Mn, 
Fe, Co, and Ni are considered and we increased the abundances of the single elements 
individually and for some combinations in the static calculations. The best-fit model 
contains about 1000 x solar nickel and about 10-100 x solar iron (Table 2). The other 
iron-group elements have almost no influence on the overall flux shape. 

Table 2. Parameters of EC 1 1481-2303. Teff, logg and the abundances for H-O 
are adopted from Rauch et al. (2010). IG denotes Ca - Mn and Co. 

Jeff/K logg [H] [Hi] [c] [N] [o] [iG] m \mr 

55 000 5.8 0.127 -1.405 -0.875 -0.875 -0.875 1-100 10-100 1000 



The stratified model yielded abundances for the iron-group element that are all su- 
persolar. This confirms the abundance results from the static model. Both the stratified 
and the static models can reproduce most parts of the observed UV flux shape. 

If the absorption lines of the lUE and FUSE observations are considered, both 
models are unable to reproduce these lines. From the fit to the Hor line, a rotational 
velocity of 30 km/ sec can be determined. Applying this to the models in the UV range 
improves the fit and some parts of the observation are matched. 

The interstellar fine absorption spectrum was modeled with the program OWENS 
and some of the strongest lines are of interstellar origin. But still, the atmosphere 
seems to be dominated by iron-group elements. To model those lines, atomic data must 
be rehable. At the moment no such data are available for temperatures in the range 
of hot sdOB stars {Tes > 40 000 K). Most lines are theoretically calculated and have 
uncertainties in the wavelength range as well as in the f- values. We suggest that these 
uncertainties are the reason for the discrepancy between the models and the observation. 
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